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SUMMARY 


The theory of sound propagation over randomly irregular, 
nominally plane terrain of finite impedance is discussed. The 
analysis is an extension of the theory of coherent scatter 
originally proposed by Biot [28] for an irregular rigid surface. 
It combines Biot's approach, wherein the surface irregularities 
are modeled by a homogeneous distribution of hemispherical 
"bosses," with more conventional analyses in which the ground is 
modeled as a smooth plane of finite impedance. At sufficiently 
low frequencies it is shown that the interaction of the surface 
irregularities with the nearfield of a ground based source leads 
to the production of "surface waves," which are effective in 
penetrating the ground shadow zone predicted for a smooth surface 
of the same impedance. 
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LIST OP PRINCIPAL SYMBOLS 


a hemispherical "boss" radius 
c speed of sound 

f frequency (Hz) 

F excess attenuation function, equation (3.9) 
G Green's function 

h see Figure 1 

k Q to/c 

% minimum distance between neighboring bosses 
n iko/Z 

n* see (2.28), (2.31) 

N mean number of bosses per unit ground area 
w numerical distance, equation (3.7) 

^m ( x mi ' 0, x m ^), center of mth boss 

Z R + iX, ground impedance 
k flow resistivity of ground 

v see equation (2.29) 

K elevation of ground from x 2 = 0 

o fractional ground area covered by bosses 
p source-observer distance along ground 
<f> velocity potential 

to = 2irf radian frequency 
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1. INTRODUCTION 


A plane sound wave is incident on a flat surface of finite 
impedance. As the angle of incidence approaches grazing the 
amplitude of the reflected wave tends to that of the incident 
wave, but their phases ultimately differ by 180° [1], Thus a 
plane acoustic wave cannot propagate parallel to the surface, 
and, in particular, this implies that most of the surface would 
be in a "shadow zone" of an acoustic source located on or very 
close to the surface. The intensity of the sound in this shadow 
zone is very much less than it would be at the same distance from 
the source had the surface been an ideal, rigid reflector. The 
point source problem has been studied extensively in the 
literature (see, e.g., the reviews [2-4]), and is the acoustic 
analogue of the problem of radiation from an electromagnetic 
dipole situated above a conducting boundary [5-7]. Detailed 
analyses for an acoustic point source are given in references [8- 
17]. 


These studies indicate that the shadow zone is penetrated by 
a "ground wave" whose properties depend on the surface impedance 
characteristics. Under certain conditions the ground wave 
includes a trapped surface wave which propagates subsonically in 
the acoustic medium just above the surface. The amplitude of 
this wave decays inversely as the square root of the propagation 
distance r, say, due to cylindrical spreading from the source, 
and exponentially with r due to absorption by the surface. If 
the absorption is sufficiently small, however, there can be a 
considerable interval in r within the shadow zone wherein the 
acoustic field exceeds that in free space at the same distance 
from the source (i.e., the surface wave provides a temporary 
means of negative damping). 
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There is confusion and controversy in the literature 
regarding the role of the surface wave component of the ground 
wave [14,16]. This has arisen because account has not always 
been taken of changes in the asymptotic formula for the ground 
wave function F(w) in different regions of the w-plane, w being 
the complex numerical distance (see Sec. 3). Accordingly, 
experimental studies often involve conflicting interpretations of 
the measured sound pressure levels (see e.g., refs [2,3,18-22]). 

The surface wave phenomenon has recently been discussed by 
Tolstoy [23-27] in connection with the propagation of sound over 
a randomly irregular surface. This work is a development of an 
earlier treatment of the problem due to Biot [28-30], involving 
propagation over a rigid surface formed by a homogeneous 
distribution of hemispherical bosses on a plane. The surface 
diffracted sound comprises a coherent component (obtained by 
averaging with respect to an ensemble of realizations of the 
distribution of bosses) and an incoherent component (having zero 
mean). When the acoustic wavelength is large relative to the 
spacing of the bosses, the incoherent field is small, and may be 
neglected. In these circumstances the coherent radiation is 
found to satisfy an impedance condition on a mean, smoothed 
plane, an important property of which is that it leads to the 
prediction of a coherent surface wave. Since the composite 
irregular surface is rigid, however, this wave is not attenuated 
by absorption at the boundary, although it must ultimately decay 
due to scattering into incoherent waves. At low frequencies the 
influence of such scattering is important only over very large 
propagation distances. In consequence there is effectively no 
surface shadow zone; on the contrary the surface pressure tends 
to be enhanced. Tolstoy [24] has demonstrated theoretically that 
the surface wave can also eliminate shadowing caused by a 
vertically decreasing sound speed profile (which tends to refract 
waves away from the surface). 
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Experiments performed by Medwin et al [31], using an 
acoustically rigid irregular surface, confirm the presence of the 
surface wave, and are in excellent detailed agreement with the 
coherent wave theory. It is possible that this effect will also 
be significant for propagation over homogeneous irregular ground 
cover of finite impedance, and that there may exist ranges in 
frequency and propagation distance within which the surface wave 
will provide a degree of negative damping. 

In this paper the Biot/Tolstoy theory is extended to examine 
the interaction of sound with an irregular surface of finite 
normal impedance. Attention is given, in particular, to the 
problem of long range propagation over irregular terrain when 
both the source and observer are on or close to the ground, and 
when, in addition, the relevant acoustic wavelengths are large 
relative to the scale of the surface irregularities. The latter 
approximation is necessary to make the analysis tractable, but 
may otherwise be justified on the basis that sound whose 
wavelength is comparable to, or smaller than, the scale of the 
surface irregularities, will be rapidly dissipated by surface 
scattering. Indeed, according to Piercy et al [2], for broadband 
noise sources and propagation distances of a mile or more, the 
high frequencies are rapidly attenuated by atmospheric 
absorption, midfrequencies by ground shadowing, and the main 
contribution to the measured sound pressure level is from surface 
waves of frequencies f < 200 Hz. 

In Section 2 of this paper the Biot/Tolstoy boundary 
condition for the coherent field is determined for an irregular, 
locally reacting surface. The analysis is restricted to one-and 
two-dimensional surface roughness elements modeled respectively 
by a distribution of parallel cylinders of semi-circular cross- 
section, and hemispherical bosses. There is, of course, no 
necessity to limit consideration to such specific roughness 
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shapes, but more general formulations of the theory (c.f. [23]) 

are unlikely to alter the principal conclusions of this paper, 
provided the distribution of surface irregularities is 
sufficiently homogeneous. The boundary condition is used to 
study propagation from a ground-sited point source to an observer 
also positioned on the ground (Sec. 3), making use of empirical 
estimates of the ground normal impedance. In particular the 
influence of surface roughness on the ground shadow zone is 
examined . 
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2. COHERENT ACOUSTIC BOUNDARY CONDITION AT A ROUGH SURFACE OF 


FINITE NORMAL IMPEDANCE 

Consider the propagation of sound in the region 
x 2 > C ( x Q ) ( a= 1,3) of a rectangular coordinate system 
x 1 ,x 2 rX 3 , where x 2 = c(x a ) the tough surface formed by a 
distribution of hemispherical bosses (N per unit area) of radius 
a over the plane x 2 = 0 (see Figure 1). The length i ~ 1//N is 
taken to characterize the spacing of the bosses, and it is 
assumed that £>>a. The sound has time harmonic dependence e lut 
(which is suppressed in the following), and satisfies the 
Helmholtz equation 

(72 + kZ)* = 0 , x 2 >c(x q ), (2.1) 

<j> being the velocity potential, k Q = co/c >0, and c is the speed 
of sound. Greek suffixes, a, 3, etc, range over the 1, 3 
directions only. 

On the irregular surface the potential satisfies 

— + n<(> = 0 , x 2 = C(x^) ' (2.2) 

n 

where x is a local coordinate normal to the surface and directed 
n 

into the acoustic medium, and where 



is the specific normal impedance , assumed to depend on frequency 
alone. Attention is confined to sound of sufficiently low 
frequency that 
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(2.5) 


ITT ~ k o a << 1 * 

The magnitude of Z in terrestrial applications is discussed in 
Section 3. 

Introduce a control surface S : x 2 = h (h> a), just above 
the irregular boundary, and consider the diffraction of an 
incident wave denoted by 4 ^ In x 2 > h the total field can be 
expressed in the form 

4(x) = ♦ I (x) + <J> R (x) + | -|y7 ( ^ } g( £'£) d yi dy 3 , (2.6) 

the surface integral being taken over S. In (2.6) <f> R (x) is the 
wave that would be specularly reflected from a rigid plane 
coinciding with S, and the Green's function G(x,^) has outgoing 
wave behavior in x 2 > h, satisfies 9G/3y 2 = 0 on S, and is given 
by 


G (x,y) 


-1 

TiF 


( ik 0 !s“£l 
) e 
[Tx^T 


+ 


ik 0 \*-Z\ 

e 



(2.7) 


where £ = (y lf -y 2 + 2h, y 3 ). 

We now derive the leading order approximation to the 
coherent (average) form of (2.6). 

The integration in (2.6) may be expressed as the sum of 

integrals over a sequence of non-overlapping area elements S m of 

S. Each S m covers precisely one hemispherical boss, and has a 

diameter of the same order as the separation distance 1 . 

Let x = (x , 0, x ) denote the position of the center of the 
m| m3 

boss in S m . in the neighborhood of x m the coherent field 
4>(x), say, can be expressed in the form 
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+ <£> - * ( im> + (x j - x mj ) 4xJ ( *. 

3 2 $ 


<?»> 


+ o(x, “X . ) (X, -X . ) 


2 i mi j mj' ax^x.'-m 


(x ) + . . . 


( 2 . 8 ) 


(no summation with respect to the repeated suffix m). It follows 
from (2.2) - (2.5) that 


3 <(> 

3x 


a 



3<j) 

3x 2 



(2.9) 


and that the final term on the right of (2.8) ~0(k Q £) 2 4> , at 
most . 


This ordering suggests the following representation 


of <J> ( x ) in S m; 


<t>(x) = <f>(x m ) + <f>i(x) + 4> 2 ( 5.) + •••' 


( 2 . 10 ) 


where 


♦i ( &> = (x a + 2|x a -x 

• ~ ~m a 


( 2 . 11 ) 


~ 0 ( k Q £ ) <J) 


and <(>2 contains all remaining terms ~0(k Q £) 2 4>. Observe that, 
correct to 0(k Q £) 2 , this representation implies that 


V 2 <j> ! = 0 ; V 2 4> 2 = -kj?(x m ) , 


( 2 . 12 ) 


and that 3<j>i/9x n = 0 on the surface S mQ comprising the portion of 
the plane x 2 = 0 and hemispherical boss in S m (see Figure 1). 
Hence, to the same approximation, (2.2) becomes in S m : 
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(2.13) 


‘ ” n ^ ( ~ma ) ' x 2 = « (x „> • 

Next, for an observer at x at a large distance 
from £ in S (2.7) becomes 

G(x, t ) = Gfx.xJ + <y o -x mo ,)-||— (x.x^) + .... (2.14) 

ma 

where each term in this expansion is 0(k Q £) relative to the 
preceeding one. 

The expansions (2.10), (2.14) are now used to determine the 

principal contribution from S m to the integral in (2.6). 

Consider first that arising from 4>i(x). Since £>>a, one finds, 

correct to the retention of terms of order (k a) 2 , 

o 


t ly7 ( £ )G( £'£ )dy i dy 3 




~ i (y cT x m<x ) ly7 ( X ) 13T ( S'*J dy l dy 3 

hi nnr» 


■’ a3 

m a ma 


(2.15) 


To the same order of approximation, the net coherent contribution 
to the surface integral in (2.6) from all terms of this type can 
be set in the form 

3 1 3 2T 

Z V 177 G{ Z'V d y i d y 3 ■ 4> oa 1777 G ^,JC)dy idy 3 , (2.16) 


where 


= Nua 2 


(2.17) 
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is the mean fractional surface area of the plane x 2 = 0 covered by 
the hemispheres. 


Similarly, making use of the divergence theorem, one obtains 


9^2 


942 


4 iyf ( £ )G( £'£ } d y idy 3 - $ -§y- (x> G <£'Sj d yi d y3 


m 


m 


G (X'*J < $ V 2 4> 2 d 3 iC + $ V $2 * d £ 

n ) V C 

l m m 


+ ^ Vif> 2 * d £ 

s' 

mo 


(2.18) 


where V m is the volumetric region bounded by S m , S mo (see Figure 
1) and laterally by the surface C m , and dz is the surface element 
directed into V m . 

The individual coherent components of the integrals in the 
brace brackets of (2.18) are evaluated as follows: 


$ v 2 «f> 2 d 3 ^ 

V 

m 

<$> v<j> 2 . dz 

c 

m 


" - k o d yi d Y3 • 

S 

m 

- -i^ ( s m )h $ d yi d v3 ? 

a S 


(2.19) 

( 2 . 20 ) 


2 • dZ « -n(l + a) $(x ) 
S 

mo 

In obtaining these results use has 
(2.13), and it has been noted that 


4 dyjdy 3 . (2.21) 

5 

m 

been made of (2.12) and 
on C m , <j> 2 may be replaced by 


" *2 -|x 2 ( s m > + 

B a 


( 2 . 22 ) 
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correct to an error (due to the field scattered by the boss in 
S m ) which is 0[(a/£,) 2 ] smaller. 

Summing these results over all area elements S m , as 
previously, combining with (2.16) and substituting into (2.6) we 
obtain the coherent field representation: 

♦ (x) = <J> (x) + 4 > d (x) +<£ |ga— 1 2 | -h-^I 

1 ~ R ~ s < 

-(h-|2a)k2 $ -n(l+a) 4> }G(x,y)dy jdy 3 (2.23) 


In other words, for the coherent component, the effective 
value of 3^/9x 2 on x 2 = h is given by 


3 cj> 
9 x 2 


- ( h-oa ) 


9 2 <j> 

9x z 
a 




+ 


n ( 1 + a) 



(2.24) 


But <j> satisfies the Helmholtz equation (2.1), so that this result 
is equivalent to 


9 <t> 
9X 2 


( h- oa ) 




n( 1+a) 



(x 2 = h) . 


(2.25) 


The first term on the right of (2.25) may be discarded, 
since it accounts for a small uniform phase shift in the 
reflected field which is of no practical significance. To see 
this, note that, if the wall were rigid and perfectly flat 
(n=0, a=0), (2.25) would reduce to 

Ixl ~ h Txf ' 0 on x 2- h ' (2.26) 

which, correct to a relative error ~ 0(k Q h) 2 , is the same as the 
exact rigid surface condition 
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t£ 2 “ 0 ° n X2_ 0 * 

In the presence of surface irregularities (2.25) indicates 
that the effective location of the smoothed plane for application 
of the boundary condition is x 2 = oa. This is very much smaller 
than the acoustic wavelength of interest, and henceforth the 
condition will be applied on x 2 = 0, in which case (2.25) reduces 
to 

= -ik Q n*$ , x 2 = 0 (2.27) 


where 


n* 


l+o 

z 


- ivk o a 


(2.28) 


and v = a/3. 


In obtaining this result it has been implicitly assumed that 
the hemispherical bosses are sufficiently far apart that mutual 
interactions via the incoherent field are negligible, an 
approximation which requires a<<l. For finite a the analysis may 
be modified (along the lines described in Biot's [28] original 
work, or as described in reference [32]) to deduce that (2.27), 
(2.28) remain valid provided that the definition of v is modified 
as follows 


v = oj 


(1 + 71' 


2 

3 


(2.29) 


where £ is the minimum distance between the centers of the bosses 


on x 2 = 0. 
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One Dimensional Surface Roughness 


For an impedance surface which is irregular in, say, 
the xp direction alone, the roughness elements may be modeled by 
a surface distribution of cylinders (N per unit length in 
the xp direction) of semi-circular cross-section and radius a. 
The long wavelength, coherent field boundary condition in this 
case becomes 


8 4 > 

3x 


= -ik n*$ - ^ a 
2 o Y 4 


3 2 4 > 
9 ^ ' 


X 2 — 0 , 


(2.30) 


where o = 2Na is the mean fractional area of the plane x 2 = 0 
covered by the cylinders, and 


n * = 1 + ( rc-1 ) a 

Z 


i it 
4 


ok a 
o 


(2.31) 


These formulae are valid for a<<l only. 
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3. LONG RANGE PROPAGATION FROM A POINT SOURCE CLOSE TO THE 
GROUND: TWO-DIMENSIONAL SURFACE ROUGHNESS. 

The formulae derived above are now applied to study the 
propagation of low frequency sound over randomly irregular 
terrain. Both the source and observer are assumed to be located 
at ground level (or close to he ground, since the wavelength is 
large). This restriction avoids many of the complications which 
tend to obscure more general treatments of propagation over an 
impedance surface. 

Consider the outgoing wave solution (for the coherent field) 
of 

(V 2 = k § ) <J> = i(x 1 )i(x 2 -e)i(x 3 ), e>0, (3.1) 

where <j> satisfied (2.27) on x 2 = 0, and the overbar on the 
coherent field is henceforth suppressed. At ground level, 
x 2 = Of and as e> + 0, the following Fourier integral 
representation is readily derived with the aid of the boundary 
conditions 

„ KJ 0 (Kp)dK 

0 — / (3.2) 

{/K^-kg -ik 0 n* } 

where J 0 is the zeroth order Bessel function, and p = 

✓x? + x§ . A branch cut for /K^-k § in the half-plane Re K > 0 
extends along the real axis from K = k 0 to + «. The path of 
integration runs along the lower side of the cut ( ImK = -0) where 
/K 2 -k g is real and positive. This configuration ensures that 
Im /K 2 -k § < 0 in the cut half-plane Re K > 0, which, for an 
absorbant, irregular surface (R > 0 in (2.4)), implies that the 
simple pole of the integrand at 


<P(P) = or / 
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K = k 0 ( 1 - n* 2 ) 1 ^ 


- k 0 (l 



( n* 2 <<1 ) 


(3.3) 


does not lie in the cut half-plane. 


The decomposition Jg(Kp) = (H 0 (Kp) + H 0 (Kp))/2 in terms 

/ 1 t ( 2 ) 

of Hankel functions Hg , Hg permits (3.2) to be replaced by 
the sum of three integrals. Two of these, the first involving 
Hg ^(kp) and the second Hg ^(Kp), are taken respectively along 
the positive/negative imaginary axes, and sum to zero ([33]), 

Sec. 9.6.4). The third integral is over the branch cut and 
reduces to 


<J>( p) 


-1 
2ir p 




(kgp/1 + X 2 )dX 
X 2 + n* 2 


(3.4) 


([34], page 736). In this expression the first term in the curly 
brackets is the potential that would be observed on a perfectly 
flat, rigid surface. For an observer in the acoustic far field 
(kgp >> 1), the value of the remaining integral may be 
approximated by use of the asymptotic formula 


H 


(!) 

0 


(Kgp/TTT 2 ) 



e i{k 0 py'l + X 2 - tt/4 } 

(l + x 2 ) 1/4 


(3.5) 


According to the principle of stationary phase the dominant 
contribution to the integral is from the region X << 1, and may 
be obtained by expanding the integrand about X = 0. In doing 
this, however, note that the pole of the integrand can lie 
arbitrarily close to X = 0 when the specific impedance 2 is 
large. To account for this the full effect of the term X 2 + n* 2 
in the denominator must be retained, leading to 


4> ( p ) 


w 


-e 


ikgP 
2ir p 


[1 


2we 

/ it 


i tt/4 


ix 2 dx 

r 00 e, 

J o x z - Iw ' 


(3.6) 
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where, in accord with conventional notation, 
ik 0 pn* 2 

w = 2 (3.7) 

is the (complex) numerical distance between source and observer. 
Since n* 2 << 1, |w| can assume all values in the range 
0 < |w| < » when k 0 p is large. Finally, a series of standard 
transformations ([33], Sec. 9) enables (3.6) to be expressed in 
the form 

-e ik ° p 

* (p) = ~ 2 tt P ' F(w) (3.8) 

where F(w) is defined in terms of the complementary error 
function by 

F(w) = 1 + i /ttw e w erfc(-i /w~ ). (3.9) 

In evaluating this expression we take 

/w = Iwl^ 2 exp{ ^ arg w}, where < arg w <3tt/2. 

This representation of the far field is consistent with 
previous studies of propagation over a plane impedance surface 
(c.f. references cited in the Introduction). The function F(w) 
(hereafter designated the excess attenuation function ) accounts 
for the influence of the irregular impedance surface on 

propagation relative to the field <J> 0 (p) = -e lk ° p /2Trp obtaining 
for a smooth rigid plane. F(0) = 1, and when |w| >> 1 , 

F(w) » 2i /ttw e w + < arg w < ~ ; 

i 3 

" 2w + * * ' 7 < arg w < T~ * ( 3 . 10a , b ) 
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The first term on the right of (3.10a) represents the so-called 
surface wave . This is exponentially damped if the boundary 
resistence R / 0, although there may exist intervals in |w| 
wherein |F| >1. In the extreme case in which R = 0 


w 


-ik 0 P 

2 


r 1 - * ~ + vk 0 a) 2 


(3.11) 


Since it may always be assumed that k 0 has a small positive 
imaginary part (corresponding to the existence of viscous 
dissipation in the acoustic medium), it follows that 
arg w = — tt/2 + 0, and therefore that, for |w| >> 1, 


<!>( p) 


0 


1 + 


2tr p 


+ vk 0 a 


e i {k 0 p ( 1 


+ 



2 ) - */4} + 


0 ( 1/wp ) 


(3.12) 


Thus, the field at ground level is dominated by the surface wave, 

1 2 

which propagates subsonically at the phase velocity c/{l + -j|n*| }, 

and decays like l//p 

For an irregular, rigid surface (X = «° ) these results are 

equivalent to those obtained by Tolstoy [23], and are in accord 

with experimental observations of Medwin et al [31]. These 

experiments confirm the predominance of the surface wave in the 

far field, and in particular that the amplitude decays like l//p 

3 /2 

and is proportional to k 0 (for X = »). They also establish 
the validity of predictions of the coherent wave theory 
for k 0 £ < 1. 
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3.1 Propagation Over Randomly Irregular Terrain 

The relative importance of irregularities in ground 
elevation on the propagation of low frequency sound will now be 
examined. The dependence of the specific normal impedance 
^ R + ix on the frequency F = co/2tt Hz is approximated by 

R = 1 + 0.051(</f) 3 / 4 

(3.13) 

X = 0.077( </f ) * 73 

< being the flow resistivity of the ground. A reasonable fit to 
experimental data given in [2,3] for grass covered ground is 
obtained when k = 230,000 MKS units. Observe that )z| — l < 0.1 
for F < 450 Hz, and does not excees 0.175 until f > 1000 Hz. 

This may be taken to justify conditions imposed above in (2.5). 
Evidently, for such frequencies the contributions to n* (Equation 
(2.28)) from Z and the surface irregularities can be of the same 
order of magnitude. 

The variation of excess attenuation, 201 og i 0 | F | dB with 
absolute numerical distance |w| is illustrated in Figures 2(a) - 
2(d) for surface irregularities of different sizes and 
concentrations for the two frequencies f = 10, 100 Hz (for which 
I Z | 1 = 0.0065, 0.035, respectively). The case a = 0.9096 

corresponds to tightly packed hemispherical bosses. The dashed 
curves denote the excess attenuation ("ground shadowing") for 
flat ground (a = 0), where ultimately <J> ~ Z 2 /k 0 p 2 as p -► «. For 
flat, rigid ground, F = 1, and this characterizes the "near 
field" behavior of the curves in the Figures for small values of 
|w|. When k 0 a is large the surface irregularities have a 
significant effect on propagation, the peaks in the curves 
corresponding to the local predominance of the damped surface 
wave term (first term on the right of (3.10a)). In such regions 
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the onset of the ground shadow zone is delayed by the surface 
wave to such an extent that the field exceeds that obtaining for 
a smooth, rigid plane. At larger distances the surface wave is 
absorbed by the ground. 


3.2 One-Dimensional Surface Irregularities 


In this case formula (3.9) for the excess attenuation F(w) 
remains valid provided the numerical distance w is defined by 


w 


ik 0 P 

~2 


{ 


( 1 + ( TT 

Z 


l)o) 


i tt ok pa 

3 


cos 2 0 


r 


(3.14) 


for o < 1. Here 0 is the angle between the direction of 
propagation parallel to the ground and the x x axis (the surface 
corrugations being parallel to the x 3 -direction ) . The variation 
of the attenuation with |w| for different propagation directions 
6 is illustrated in Figs. 3(a), (b) for o = 0.2, k(ja = *1 and for 
f = 10,100 Hz respectively. At 6 = 90° the propagation charac- 
teristics reduce to those for smooth ground of specific normal 
impedance Z/(l + (it - l)o). 
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4. CONCLUSION 


The theory of low frequency sound propagation over a 
statistically homogeneous irregular surface of finite impedance 
has been developed by extension of the Biot/Tolstoy [23-30] 
theory of coherent scattering by an array of rigid bosses 
distributed over a plane. In terrestrial applications, the 
approximations involved are expected to be adequate at 
frequencies less than about 1000 Hz. The presence of the surface 
irregularities is responsible for the generation by a ground- 
based acoustic source of a surface wave which can be effective in 
penetrating the ground shadow zone predicted theoretically when 
the ground is modeled as a smooth, plane impedance surface. 
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FIGURE 2a. Variation of excess attenuation with absolute numerical 

distance for irregular terrain of finite impedance modelled 
by a distribution of hemispherical bosses of radius a. 

(a) a = 0,2, f = 10 Hz; 

: excess attenuation for flat ground (a = 0) at 

the same frequency. 
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FIGURE 2d .a = 0.9096, f = 100 Hz 





FIGURE 3a. Variation of excess attenuation with absolute numerical 

distance for irregular terrain of finite impedance modelled 
by a distribution of parallel cylinders of semi-circular 
cross-section. 0 is the angle between the propagation 
direction and the x^axis (normal to the axes of the 
cylinders); a = 0.2, k 0 a = 0.1: 


( a ) f = 10 Hz ; 




FIGURE 3 
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